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The almost periodic eigenvalue problem described by the 
Harper equation is connected to other classes of quasiperiodic 
behaviour: the dissipative dynamics on critical invariant tori 
and quasiperiodically driven maps. Firstly, the strong cou- 
pling limit of the supercritical Harper equation and the strong 
dissipation limit of the critical standard map play equivalent 
role in the renormalization analysis of the self-similar fluctu- 
ations of localized eigenfunctions and the universal slope of 
the projected map on the invariant circle. Secondly, we use 
a simple transformation to relate the Harper equation to a 
quasiperiodically forced one-dimensional map. In this case, 
the localized eigenstates of the supercritical Harper equation 
correspond to strange but nonchaotic attractors of the driven 
map. Furthermore, the existence of localization in the eigen- 
value problem is associated with the appearance of homoclinic 
points in the corresponding map. 



I. INTRODUCTION 

Systems with two competing periodicities have been 
in the forefront of research in nonlinear dynamics as well 
as in condensed matter physics. Perhaps the most well- 
known paradigm in the study of quasiperiodic dynamics 
in autonomous dissipative and Hamiltonian systems is 
the standard map with and without dissipation A 
related problem in the condensed matter physics is the 
Schrodinger equation with quasiperiodic potential [Q . A 
simple discrete prototype in this field is the Harper equa- 
tion H or the almost Mathieu equation Q . Invariant cir- 
cles in the standard map correspond to extended states 
of the Harper equation while the breakup of an invariant 
circle is analogous to the localization transition in the 
Harper equation. Both transitions have been studied by 
the renormalization methods . 

We have applied the decimation methodology, devel- 
oped earlier to show the existence of " slope" universality 
in the dissipative, critical standard map p[, to a variety 
of eigenvalue problems which can be written in the linear- 
difference form JTo| ] . Very recently, the renormalization 
analysis was extended to the Harper equation beyond 
criticality These studies have shown the existence 

of a unique strong coupling fixed point which describes 
the fluctuations of localized eigenstates in the whole of 
supercritical region. In the first part of this paper, we 
show that this strong coupling fixed point of the Harper 
equation is related to the strong dissipation limit of the 



critical standard map. The renormalization analysis for 
the critical standard map is rewritten in the form which 
makes the correspondence clear. 

The second part of this paper is to relate the results 
of the supercritical Harper equation to the dynamics of 
quasiperiodically forced dissipative maps. In parallel to 
the study of quasiperiodic dynamics in autonomous dissi- 
pative and Hamiltonian systems there have been investi- 
gations on the dissipative dynamics of quasiperiodically 
forced systems. Although the presence of such forcing 
causes some usual orbits like the periodic ones to be ab- 
sent in these systems, there are other interesting features 
which have not been found in autonomous systems. Per- 
haps the most striking phenomenon is the appearance of 
strange nonchaotic attractors (SNA) [Q. These are ge- 
ometrically strange attracting sets for which all the Lya- 
punov exponents (except the one related to the forcing) 
are negative. Inspite of a great deal of research in sys- 
tems exhibiting SNA |^3|,Q, there are still many open 
questions: the condition for the existence of SNA as well 
as a general renormalization theory on their appearance 
are not fully understood. 

Bondeson et al. ]l5| ] pointed out an interesting connec- 
tion between a quasiperiodically forced oscillator and the 
Schrodinger equation with quasiperiodic potential. Using 
the substitution introduced by Priifer |l6| already sev- 
enty years ago, the two classes of systems can be mapped 
onto one another the energy eigenvalue of the Schrodinger 
equation appearing as a parameter for the oscillator. The 
transition from an invariant circle to SNA in the classi- 
cal problem is related to the transition from extended to 
localized states in the quantum problem. 

In this paper, we map the discrete Schrodinger equa- 
tion, namely the Harper equation, to a quasiperiodically 
driven one-dimensional map so that the eigenfunction 
corresponds to an attractor of the map. The simplic- 
ity of the discrete system compared to the continuum 
model offers the possibility of obtaining further insight 
into the existence and characterization of SNA. By fo- 
cusing on the strong coupling limit, we are able to show 
the existence of SNA in quasiperiodically driven dissi- 
pative maps. Furthermore, this approach also helps in 
understanding some of the open questions related to the 
significance of the critical phase in the Harper equation. 
HjyJ Our focus is somewhat similar to the one of Bon- 
deson et al. [fl5| in that we are interested in various tran- 
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sitions and their effects on the renormalization dynamics 
or the dynamics of the forced map. 

In Section II we briefly review the renormalization ap- 
proach as applied to the Harper equation and the dis- 
sipative standard map. In Section III, we show how 
these two problems are related. In Section IV, we dis- 
cuss the dynamics of the one-dimensional map associated 
with quasiperiodic eigenvalue problems and for the first 
time show explicitly the connection between localization 
and strange nonchaotic attractors conjectured earlier by 
Bondeson et al. [jl5|]. In Section V, we discuss our re- 
sults from a general perspective. As a comparison with 
previous two quasiperiodic systems, we describe the re- 
sults with completely different properties for the phonon 
equation of the Frenkel-Kontorova model although this 
eigenvalue problem is formally very close to the Harper 
equation. 

II. RENORMALIZATION APPROACH TO 
QUASIPERIODIC SYSTEMS: THE HARPER 
EQUATION AND THE DISSIPATIVE 
STANDARD MAP 

We have developed a decimation approach to study 
quasiperiodic systems with two incommensurate frequen- 
cies which can be written in the following linear-difference 
form: §£l| 

a(i)ip i+ i + b(i)ipi-i + c(i)ipi = 0. (1) 

Here, a, b, and c are real or complex functions of ia 
where er is an irrational number equal to the ratio of 
two frequencies. Although the decimation method can 
be implemented for any irrational a, here we consider 
the simplest case where a is given by the inverse golden 
ratio a = (y/5— 1)/2. Our formulation can be generalized 
also to the case where ipi is a multi-component vector and 
a, b, and c are matrices. 

In the decimation scheme, it is appropriate to decimate 
out all sites except those labelled by the Fibonacci num- 
bers F n (which are the best rational approximants of the 
golden ratio). At the n th decimation level (n = 2,3, ...), 
the linear-difference equation is expressed in the form 

f n (i)ip(i + F n+1 ) = ip(i + F n ) + e n (i)ip(i). (2) 

The additive property F„ +1 = F n +F n _i of the Fibonacci 
numbers provides exact recursion relations for the deci- 
mation functions e n and f n : 

e " +l(i) = -TTaJJI) (3) 

, f n -i(i + F n )f n (i + F n ) 

fn+l{l) = YTAfjl) (4) 

A = e n -i(i + F n ) + f n -i{i + F n )e n (i + F n ). 



The decimation functions e„ and /„ define a renor- 
malization flow which converges asymptotically on an at- 
tractor. In the case when the attractor is a limit cycle 
of period p, it reflects the translational invariance of the 
function ijj in the Fibonacci space and also captures the 
self-similarity underlying it at all scales. The decimation 
functions e n and f n determine the universal scaling ratios 

Q = lim \^(F pn+j )/^(0)\; j = 0,...,p-l. (5) 



The two quasiperiodic systems that we discuss here are 
the Harper equation 

V'i+i + tpi-i + 2A cos[27r(icr + 4>)]ipi = Eip., (6) 

and the dissipative standard map 

k 

x i+ i = Xi + + b(xi - Xi-i) - — sin(27raA (7) 

In the Harper equation, E is the eigenvalue correspond- 
ing to the eigenfunction For irrational cr, the equa- 
tion describes a particle in a periodic potential incom- 
mensurate with the periodicity of the lattice. On the 
other hand, the standard map with < b < 1, and for 
certain values of £1, describes quasiperiodic dynamics on 
an invariant circle. This equation can be written in the 
linear-difference form by considering the equation for the 
derivatives & = dxi/dxo describing the slope of the i 
times iterated reduced circle map Xk+i = h(xk) at xq\ 

6+i + -[l + b-k cos(27rs i )]6 - 0. (8) 

In ref. || , a somewhat similar formulation of the dec- 
imation equation (||) and the recursion relations was 
used to derive the following result along the critical line 
k = fc c (6, f2) for the breakup of the golden invariant cir- 
cle in the standard map: £,F n b~ Fn — * Cs as n — * oo, 
where (s — 0.435625... and the slope was calculated at 
the point x Q which lead to the smallest value. In fact, the 
value of £ was calculated in the limit 6^0 and its uni- 
versality, termed as the "slope universality" , for finite b 
was not understood although suggested by the numerical 
results. This result, which has not been yet rigorously 
proven, would have an interesting mathematical implica- 
tion: under renormalization, the critical renormalization 
fixed point would attract not only circle maps with zero 
slope but also those with finite slope. The latter would 
correspond to critical reduced circle maps of higher di- 
mensional systems with finite Jacobians. They would 
give a new dimension for the stable manifold of the crit- 
ical renormalization fixed point. 

Later the same decimation procedure was applied 
to the critical Harper equation with A = 1 for fixed E 
which was usually chosen either at the band center or at 
the band edges. Here the p universal numbers £ described 
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the amplitude of the wave function at points which were 
spaced p (or any multiple of p) Fibonacci sites apart from 
the central peak. More recently the modified deci- 
mation procedure explained in above was used to describe 
universal features in the Harper equation beyond critical- 
ly. Q Here, the fluctuations of localized eigenfunctions 
were shown to be universal. 

In the next section, we discuss both problems using 
this modified decimation approach and show that the 
supercritical Harper equation with A > 1 and the critical 
dissipative standard map with b < 1 are closely related. 
In particular, this analysis clarifies the universality of (s 
for the full range of the dissipation parameter b. 

III. ANALOGY BETWEEN THE 
SUPERCRITICAL HARPER EQUATION AND 
THE CRITICAL DISSIPATIVE STANDARD MAP 

In this section, we explain the common principle be- 
hind the slope universality of the critical dissipative stan- 
dard map H and the supercritical universality of the 
Harper equation JTT[ . In both cases, the universality can 
be generalized also to other than the Fibonacci sites but 
for the present purpose it suffices to consider this re- 
stricted universality. 

The Harper equation in the supercritical regime (A > 
1) is characterized by the positive inverse localization 
length 7 = ln(X) for the wave function ip* |l7j]. The 
fluctuations of the localized eigenfunction are obtained 
by factorizing the exponentially decaying part out. The 
wave function ipi is written as 

Vi = e-^V (9) 

The resulting equation for the fluctuations is given by 
the linear-difference form 

\ni+i + A?ft_i + 2A cos[27r(icr + <p)]r]i = Er^. (10) 
A 

A correspondent equation for the standard map is ob- 
tained by factorizing out the exponentially decaying ef- 
fective Jacobian from the slope. In other words, with 
£i = b l rji we obtain 

brii+i + Vi-i - [1 + b - kcos(2Trxi)]f]i = 0. (11) 

It is interesting to note that the above equations have 
well-defined limits when A — > oo or b — > 0, respectively. 
Moreover, these limits resemble one another. As the lat- 
ter corresponds to the circle map, for which the existence 
of a renormalization fixed point has been known for a 
long time the analogy suggests the existence of a 

similar renormalization fixed point also for the fluctua- 
tions in the Harper equation above criticality. In fact, 
applying the decimation tools to Eq. ([ll]), the fluctua- 
tions of states with eigenvalue either at the band center 



or at the band edges were found to be self-similar charac- 
terized by periodic attractors of the renormalization flow 
PJ. In particular, for the states at the band edges with 
4> = 1/4, |?/ Fn /?7o| -> Ch = 0.172586410945... as n -> oo, 
independently of the value of A > 1. This universality 
can also be interpreted in an alternative way: the finite- 
size inverse localization length 7, defined by ipi — e~^ 1 , 
varies around the asymptotic value 7 in a universal way 
given by the equation 

7 -7 = log|^J/F„. (12) 

This is a rather intriguing result as the inverse localiza- 
tion length does depend upon A. 

In order to understand the origin of the parameter- 
independent universality in the Harper equation and the 
dissipative standard map, we write the equation for r/i in 
both cases in the form 

f 2 (i)r)(i + 2) = v(i + 1) + e 2 {i)v(i) (13) 

which gives the first step of a decimation approach of 
generating equations of the form 

f n (i)v(i + Fn+i) = V{i + F n ) + e n (i)r)(i). (14) 

The power of this method lies in the fact that in iter- 
ating the recursion relations the decimation function f n 
asymptotically "renormalizes" to zero. This is a numer- 
ical observation but also suggested by the appearance of 
the constant C — 1/A < 1 (for the Harper equation) or 
C = b < 1 (for the standard map) in front of the function 
f%. Because one can take the trivial conditions /1 = 1 
and e\ = 0, in /„ there appears a multiplying constant 
C^" -1 which tends to zero as n — > 00. In other words, 
asymptotically the recursion relation for the decimation 
function e n simplifies into 

e n+ i{i) = -e n -i(i + F n )e n (i) (15) 

which is the form in the strong coupling limit A — * 00 
of the Harper equation or the infinite dissipation limit 
b — > of the standard map. We can therefore expect 
the renormalization fixed points in these limits to de- 
scribe the scaling properties of r/(i) for the whole A or 
b range. One should also note that because /„ vanishes 
asymptotically, the scaling ratio £ is given simply by the 
asymptotic value of e n at the lattice site i = 0. 

In order to solve for the fixed points, the discrete lattice 
index i has to be replaced by a continuous variable. For 
the standard map, we can take this to be simply the 
original variable x so that the recursion relation becomes 

e n+ i(x) = -e n -i{h Fn (x))e n (x). (16) 

This equation appeared also in the previous work [Bj. 
However, because there a variant of the decimation 
method was applied directly to Eq. (pi) and not to Eq. 
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(|l l|) , one was left with another non-trivial decimation 
function for positive b. This hindered recognizing the role 
played by this equation for finite Jacobians. Because of 
the universal scaling properties of the reduced circle map 
h, e n (xo) approaches a universal limit which can be ob- 
tained in terms of the universal fixed point of the circle 
map renormalization without solving e(x) for arbitrary 
x However, one has to be careful in the analysis be- 
cause effectively it means calculating the limit of £f„ /b Fn 
taking both the numerator and the denominator to zero 
(and n to infinity). 

In the case of the Harper equation, the continuous vari- 
able for e n is obtained from the fractional part {io~} of ia, 
x = (— er) n {ier}, and the resulting fixed point equation 

e*{x) = -e*(a 2 x + a)e*{~ax) (17) 

is solved by standard expansion methods jnj. We show 
this universal function in Fig. 1 which illustrates the 
nontriviality of the strong coupling fixed point as com- 
pared to the trivial weak coupling fixed point of 
the Harper equation. As seen in the plot, the function is 
smooth and finite almost everywhere. 

It is interesting that the above fixed point equation 
( |l7| ) was independently found by Kuznetsov et al. |14j in 
studying the birth of a strange nonchaotic attractor in 
a quasiperiodically forced map. In the next section, by 
using a discrete version of the Prufer transformation, we 
relate the Harper equation directly to such a map. 

IV. LOCALIZATION AND STRANGE 
NONCHAOTIC ATTRACTORS 

The Prufer substitution that was used to map 
the continuous Schrodinger equation to overdamped, 
quasiperiodically driven oscillator has motivated us to 
consider a similar transformation also in the discrete case. 
The main idea underlying the continuous transformation 
is to define new variables p and a in terms of the wave 
function ip = pcos(o;) and its derivative ip' = psin(a) 
so that one obtains a first order equation for a where p 
does not appear and where the potential of the quantum 
problem appears as a driving term. An analogous idea 
in the case of the discrete Harper equation would be to 
set e.g. ipi-i = piCos(oti), Vi — PiSm(ai). In fact, the 
Harper equation gives trivially a very simple equation for 
Xi = cot(a-j) = V^-i/Vv 



l+ Xi -E + 2A cos[27r(icr + $)] ' K ' 

It is clear that higher iterations of this map generate 
continuous fraction expansions in terms of subsequent 
potential terms V(i) = Acos[27r(i<7 + </>)]. In the follow- 
ing, we study this one-dimensional map and verify some 
of the conjectures of Bondeson et al. Il5) in the discrete 



context. Although the transformation above is extremely 
simple, it provides a useful tool to study various eigen- 
value problems. 

Numerical study of the map ( |l8| ) shows that the ex- 
tended phase of the Harper equation corresponds to in- 
variant curves for the map. There are infinity of such 
curves depending upon the initial condition. The tran- 
sition to localization in the Harper equation corresponds 
to transition from invariant curves to an attractor. Fig. 
2 shows the numerically obtained Lyapunov exponent of 
the map as a function of the parameter E. As seen from 
the figures, the eigenenergies are singled out in this plot 
by giving rise to Lyapunov exponents larger than those 
of forbidden values of E. Furthermore, the Lyapunov ex- 
ponent varies non-smoothly as the parameter E is taken 
from allowed and to forbidden values. This in princi- 
pal, could provide a new method to calculate the allowed 
energy spectrum of the eigenvalue problem. 

As seen in Fig. 2, below criticality, the allowed en- 
ergies give rise to quasiperiodic orbits (invariant curves) 
with zero Lyapunov exponent while the states in the for- 
bidden regime are characterized by negative Lyapunov 
exponents. In the superctical regime, the Lyapunov ex- 
ponent is always negative. At criticality, where the spec- 
trum is conjectured to be singular continuous [Q, the plot 
clearly illustrates the fact that the allowed energies form 
a Cantor set of zero measure. Another interesting aspect 
of the driven map is that the inverse localization length 
of the Harper equation gives the absolute value of the 
Lyapunov exponent. Since the localization length of the 
system is exactly known, this is an example of a system 
whose Lyapunov exponent is known exactly. 

Finally, comparing Figs. 2 (a) and (b), we see a very 
interesting manifestation of the self-duality [0 of the 
Harper equation which can be written as 

j(X,E) = -ln(X)+ 7 (l/X,E/X). (19) 

Fig. 3 shows the attractor of the mapping ( |l8| ) at and 
slightly away from criticality. The self-similarity of the 
Harper wave function at critically implies that the at- 
tractor is fractal with structure at all scales. Interesting 
aspect of this attractor is the fact that it has a zero Lya- 
punov exponent. Therefore, the map exhibits a new type 
of SNA where the two nearby trajectories converge on the 
attractor, not exponentially but with a power law. In the 
supercritical regime, the attractor is associated with neg- 
ative Lyapunov exponent which appears to smooth out 
the fragmented structure of the attractor. However, as 
we will argue later in this section, the attractor is in fact 
strange throughout the supercritical regime. 

Another intriguing result emerging from this mapping 
is that in the localized regime, the critical phase <fr c of the 
potential leading to localization around the site i = can 
be found as a "homoclinic" point where the forward and 
backward attractors of the map intersect (see Fig. 4). 
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It is interesting to note that independently of the value 
of A > 1, a homoclinic point is always found at <j> c = 
1/4. In the strong coupling limit, the homoclinic points 
coincide with the divergences seen in the attractor. This 
point of view brings out the significance of the critical 
phase factor in the Harper equation which has not been 
understood in the past. Unlike the supercritical regime, 
at the critical point, the forward and backward attractors 
seem to overlay. 

In the localized regime, the attractor is strange but 
nonchaotic. The appearance of the fractal structure can 
be easily seen coming down from the strong coupling limit 
A — > oo. Here we give a rather crude analysis of this 
problem showing just the main ideas of the argument. 

Let us first rewrite the map in terms of the variables 
x and 9 assuming E = 0: 



-1 



x + 2A cos(2tt6>) 
= 9 + a mod 1 . 



(20) 
(21) 



Note that the map remains invariant under the trans- 
formation x — > —x, 9 — > 9 + 1/2 so we can expect the 
attractor to have this symmetry as well. Therefore, in 
the following we assume implicitly that 8 is considered 
mod 1/2. 

To the first order in 1/A, the attractor, whose invariant 
measure is always uniform in 9, can be written as 



x{8) 



-1 



2Acos[27r(0-cr)]' 



(22) 



However, the above form suggests a singularity appear- 
ing asymptotically as A — > oo at 8% — 1/4 + a. This 
first order singularity gives rise to higher order singulari- 
ties which are absent in the above approximation for the 
attractor. In order to see the appearance of the second 
order singularity, note that according to Eq. fl22]), \x(&)\ 
takes values in the interval [l/(2A),oo) when 8 is var- 
ied in [0, 1/2). Thus, it is possible to find 8\ such that 
x{8\ — a) + 2A cos(27r6>2) = which gives rise to a sin- 
gularity around 82 = 9\ + a. Moreover, asymptotically 
82 — > 81 + a because the singularity around 8\ becomes 
very sharp and the attractor is mostly (as far as the in- 
variant measure is concerned) close to zero away from 8\ 
as A — > 00. In other words, we have generated a new sin- 
gularity close to 82 = 1/4 + 2(7. Looking at the explicit 
(asymptotic) equation for 8%, 



4A 2 cos(27T0|) cos[27r(0| - a)} = 1, 



(23) 



shows that the second order singularity is even sharper 
than the first order one. Now the same argument can 
be repeated to generate a third order singularity from 
the second order one and so on. Each new singularity is 
sharper than the previous one. Because the measure is 



uniform in 8, the high order singularities become invis- 
ible for finite number of iterations of the map. Asymp- 
totically with A — > 00, the singularities generated in this 
way appear at m = 1/4 + mer (m = 1, 2, ...). For finite A, 
their locations move a little bit but we still expect them 
to be dense in 9. Their existence causes the attractor to 
be nowhere differentiable and to have the characteristics 
of a SNA g]. 

In summary, the above reasoning provides a very 
strong argument for the existence of SNA in quasiperi- 
odically driven maps given by Eq. (fl8|). 



V. DISCUSSION 

In this paper, we have shown an interesting relation- 
ship between the supercritical Harper equation and the 
critical dissipative standard map: the fluctuations in the 
localized eigenstates of the Harper equation are related 
to the tangent orbit of the standard map where the pa- 
rameter 1/A plays the role of the dissipation parameter 
b. In particular, the strong coupling limit of the Harper 
equation is analogous to the strong dissipation (or circle 
map) limit of the standard map. Inspite of these similar- 
ities, the two problems are quantitatively different due 
to the fact that they are characterized by different uni- 
versal numbers £. The roots for these two quantitatively 
different universality classes is tied to the fact that the 
dynamics governing the nonlinear potential in the Harper 
equation is a pure rotation while the dynamics underlying 
the critical dissipative standard map is highly nontrivial. 

The same kind of reasoning can be applied also to an- 
other related class of problems studied by us recently 
p8| , namely the phonon modes of the Frenkel-Kontorova 
model where the locations of the particles are described 
by the area-preseving standard map. Although this prob- 
lem bears a close resemblance to the Harper equation, 
the supercritical regime in the Frenkel-Kontorova model 
has very different characteristics. Below criticality, the 
potential is derived from the dynamics analytically con- 
jugate to the pure rotation in both models. However, 
beyond criticality, the Harper potential is still given by 
the pure rotation whereas the phonon potential is ob- 
tained from cantorus dynamics which cannot be mapped 
smoothly to the rotation. In the Harper equation, the 
"supercritical" regime is described by exponentially lo- 
calized eigenfunctions with universal self-similar fluctu- 
ations characterized by a unique strong coupling fixed 
point In contrast, the phonon eigenmodes defy lo- 
calization and remain critical with scaling characterized 
by a line of renormalization limit cycles. In addition, 
there exist an infinite sequence of parameter values in 
the supercritical regime where the renormalization limit 
cycle degenerates into a trivial fixed point. A very in- 
triguing characteristics of these parameter values is the 
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fact that the corresponding phonon eigenfunction is rep- 
resented by an infinite series of step-functions. 

Our studies suggest that the solutions of eigenvalue 
problems where the potential is determined by sequences 
with complicated dynamics can be extremely rich. By 
studying the Harper equation and the phonons in the 
Frenkel-Kontorova model, we may have barely scratched 
the surface of the wealth of novel phenomena exhibited 
by these types of systems. For example, in the general- 
ized versions of the Harper equations one encounters fat 
critical regimes with ergodic renormalization dynamics 
in the strong coupling limit [|lO],[ll) . 

The idea of relating an almost periodic eigenvalue 
problem to a quasiperiodically driven map has opened 
new possibilities for studying SNA. In this paper, we have 
shown that the existence of localization in the eigenvalue 
problem is associated with the appearance of homoclinic 
points in the corresponding map. Furthermore, the crit- 
ical state at the onset to localization is found to corre- 
spond to a new type of strange attractor with zero Lya- 
punov exponent. In future, we also hope to answer the 
question about the class of models whose whole super- 
critical phase is described by a unique strong coupling 
limit of the theory. In particular, it would be interesting 
to know whether quasiperiodically driven maps studied 
earlier exhibit this type of universality or not. 

The mapping from linear-difference equations to driven 
maps is of continuing interest to us. In future, we would 
like to apply this technique to a variety of other problems 
that can be written in the linear-difference form. In par- 
ticular, the one to two- hole transition in the extended 
standard map |l9| ], which corresponds to nonexponen- 
tially localized phonons [l8[ ], may provide new dynamics 
in the forced map. Fishman et al. [^0| have shown that 
the quantum dynamics of the kicked rotor can be mapped 
to the tight binding model of the above form. The inte- 
grable rotor with quasiperiodic potential has been shown 
to exhibit Harper-type universality. However, according 
to our preliminary study on the nonintegrable rotor with 
pseudorandom potential, the localized phase is not char- 
acterized by a unique strong coupling fixed point in this 
model. 

The research of IIS is supported by a grant from Na- 
tional Science Foundation DMR 093296. JAK would like 
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invitation. 
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FIG. 1. The renormalization fixed point for the fluctua- 
tions of the supercritical Harper equation obtained by solving 
the expansion of l/e(:r) upto the order x 23 by the Newton 
method. The fact that there has to be a singularity (or zero 
point for 1/e) can be seen easily from the fixed point equation. 
The horizontal line shows the value of the universal (h- 
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FIG. 2. Lyapunov exponent 7 vs the parameter E for the 
one-dimensional map derived from the Harper equation. The 
figure obtained with A = 0.5 (a) is equivalent to the figure 
with X — 2 (b) after simple scaling of 7 and E which illustrates 
the self-duality of the Harper equation. Fig. (c) shows the 
corresponding results at A = 1. The allowed values of the 
energy are characterized by an attractor with zero Lyapunov 
exponent. 
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FIG. 3. The attractor of the one-dimensional map (E = 0) 
at the critical point A = 1 (a) and at A = 1.05 (b). 



FIG. 4. The forward (lighter dots) and backward (darker 
dots) attractors of the map and the existence of homoclinic 
points for A = 2. As A approaches oo, the divergences coincide 
with the homoclinic points. 
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Fig2(a) Ketoja and Satija 




Fig2(c) Ketoja and Satija 




Fig3(b) Ketoja and Satija 




Fig4 Ketoja and Satija 




